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Abstract 

</S ■ 

■ The B-meson wavefunction has been studied with the help of the vacuum-to-meson matrix 

o : 

, element of the nonlocal operators in the heavy quark effective theory. In order to obtain the 



o 



Wandzura-Wilczek-type B-meson wavefunction, we solve the equations which are derived from the 
equation of motion of the light spectator quark in the B meson by using two different assumptions. 



. Under the condition that loq = 2A, the solutions for the B-meson wavefunction in this paper agree 
> , 

ijr^ , well with the one obtained by directly taking the heavy quark limit oo. Our results show 

(N 

CN ' that the equation of motion of the light spectator quark in the B meson can impose a strong 



constraint on the B-meson wavefunctions ^'±(c<-', z^). Based on the obtained results, we claim that 
both its distribution amplitudes (psi^) and (psi^) are important for calculating the B meson 
CLi! decays. 
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I. INTRODUCTION 



Along with the theoretical and experimental progresses, B physics is attracting more and 
more attentions. The non-perturbative light-cone (LC) wavefunction/distribtuion amplitude 
(DA) of the B meson plays an important role in making reliable predictions for exclusive B 



ummm 



meson decays. The B meson DA has been investigated in various approaches 
and is the basis of the collinear factorization shows that the B meson DA is not 

normalizable. Such feature of the B meson DA does not cause a problem in practice j^, but 
it does introduce an ambiguity in defining the B meson decay constant fs- Recently, Ref. 8| 
claims that it is the B meson wavefunction that is more relevant to the B decays and in the 

Q 

framework of the fc^r-factorization theoremlQ], they proved that the B meson wavefunction 
is renormalizable after taking into account renormalization-group evolution effects. In the 
/cT-factorization theorem, by taking into account the transverse momentum dependence (fc^- 
dependence) into the non-perturbative wavefunctions and the hard scattering, the endpoint 
singularity (an example can be found in Ref. 1^) coming from the collinear factorization can 
be cured. Theoretically, it is an important issue to study the longitudinal and transverse 
momentum dependence of the B wavefunction, since it provides a major source of uncertainty 
in the calculations of the B decays. 

Ref.^l presents an analytic solution for the B-meson wavefunction, which satisfies the 



constraints coming from the equations of motion and the heavy-quark symmetry They 
find that the "Wandzura-Wilczek-type" contribution^^ (WW approximation), which cor- 
responds to the valence quark distribution, can be determined uniquely in analytic form in 
terms of the "effective mass" (A) of the meson state, which is defined in the Heavy Quark 
Effective Theory (HQET)Ql- However, in Ref. 11 1, two extra constraints for the B-meson 



wavefunction come from the heavy quark limit, — > oo. Since the mass of b-quark is 
limited, such condition might be too strong, and we will not take such limit in our present 
calculation. In the following, we shall solve the two equations that are derived from the 
equation of motion of the light spectator quark on the basis of some physical considerations. 
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II. EQUATIONS UNDER THE WW APPROXIMATION 



In HQETfl4'], the wavefunctions "^±{1, z^) of the B meson can be defined in terms of the 
vacuum-to-meson matrix element of the nonlocal operators: 



{o|?(z)r;i,(o)|B(p)) = -I^Ti 



2t 



Here, = (0, z~, zj_), z'^ = — z^, v'^ = 1, t = v ■ z, and p'^ = Mv^ is the 4-momentum of the 
B meson with mass M. hy{x) denotes the effective 6-quark field. F is a generic Dirac matrix. 
The path-ordered gauge factors are implied in between the constituent fields. Note that in 
the above definition, the separation between the quark and the antiquark is not restricted 
on the LC {z^ = 0). 

The effective mass (A) is much smaller than the B meson mass, so the light spectator 
quark in the B meson can be treated as on mass shell with high precision. Based on the 
QCD equation of motion for the light spectator quark, we can obtain a set of equations for 
'^±{t,z'^) under the WW approximation, i.e. 

d'^-{t,z^) ^+{t,z^) -^_{t,z^) z^ d 



dt t t dz^ 

and 



[vl/+(t,z^)-vl/_(t,z^)] = 0, (2) 



d^4t,z') d^.jt^z') M/+(t,z^)-M/_(t,z^) , ^^ d^+iuj,z') _^ 
dt dt t ^ dz^ ~ 

When taking the LC limit z^ — > 0, the above two equations agree well with the ones in 



Refs 



By doing the Fourier transformation, '^±{t,z ) = J du e \E'±(u;,2; ), Eqs.(j2 



become that 



and 



Here, ujv~^ has the meaning of the LC projection of the light-antiquark momentum in the 
B meson. The exact solution of Eqs. ()4l5p will give strong constraints on the longitudinal and 
transverse momentum dependence of the B wavefunction under the WW approximation. 
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III. SOLUTIONS FOR THE B WAVEFUNCTION 



At present, one can not obtain the exact solution for the B-meson wavefunctions \I/±(c<j, z'^) 
only with Eqs. (j4l5p . Some prescriptions must be made. Under the prescription that 
z"^) have the same ^^-dependence and by taking the heavy quark limit mf, — > oo, 
Ref. [ll| obtained an analytic solution for "^±(10, z"^). In the present paper, we do not take 
such limit and try to solve Eqs. ()4l5|) by taking two other prescriptions. 

Before solving Eqs. ()4l5|) . we define two functions 4'±{uj) in the following way, 

0±H = lim ^±{uj,z^) (6) 



+0 



where the second line shows that (f)±{uj) are precisely the B meson DAs. Taking the LC 
limit {z'^ —>■ 0) in Eq.(@]), one can directly obtain a relation between the two DAs <f)+{uj) and 
of the B meson, 

.^+^,M=0, (7) 



which agrees with the one obtained in Refs. 

Furthermore, one can prove that if the B-meson wavefunctions \l/±(ct;,z^) can be con- 
structed in the following way, 

^±{uj,z^) = M^)-x{^,z^) , (8) 

then x(ci;,2;^) must be a function of the correlated variable [z'^uj{ujo — uj)Y, where cjq is the 
maximum value of u. The ansatz (jH)) shows that the two wavefunctions "$±{00, z"^) have 
the same ^^-dependence xi^,^^)- A boundary condition for xi^y^"^) can be derived from 
Eqs. fl6|8|l . i.e. lim22_,o x(c<j, 2^) = 1. From the above ansatz, we can solve Eqs. ()4p5|l exactly 
without taking the heavy quark limit rrih 00. 
Substituting Eq.® into Eq.(jD), we obtain 

x{u;,z')=auj-<P4Lu)-z')^ay), (9) 



^ The special case that x(w, z^) has nothing to do with the variable u) has been discussed in Ref.jTi] and will 
not be discussed here. 
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where ^{y) is the function of a single variable y = [uj ■ (j)-{uj) ■ z^\. The explicit form of 0-(co') 
can be derived from Eqs. ()5l9|) . which reads 

= aio^ + ao, «i = -^7^^ (10) 

where the parameter oq is to be determined. With the help of the above results, Eq.(jni) 
changes to 

with a{y) = Ay^"{y) +4^'(y) — ^(y). Eq. (fTT|) should be satisfied for all values of cu, and then 
we have to set, 

a^a{y) + (1 + a^)ay) = 0. (12) 
The non-trivial solution of the above equation is 

m = M^/y/^i), {ai < 0) (13) 

where Jo is the zero-th normal Bessel function. Here oi must be less than zero, because if 
Oi = 0, then it will lead to 0+(co') = 0; and if ai > 0, then one may obtain a self-contradictory 
result from Eq.^: ai = -i^o(0)/(4ir^(0)) = 0. 

After doing the Fourier transformation for the transverse part, \I^±(ti;, k^) = 
/ (Pzj_ exp(— ikx ■ z_|_)\l/-|-(c(j, 2;^)/(2'7r)^, the final solution for the B-meson wavefunctions can 
be written as, 

? / 1 \ , / (ki - Co'(to'o - cj)) ,^ ^, 

^±(c^, kx) = M^)^^ — (14) 

where ujq = (— ao/ai) is the maximum value of u and the normalized DAs take the form, 

2 2 
(f)+{uj) = —9{uJo - uj)uj; (l)^{uj) = —9{uJo-uj){uJo-uj). (15) 

The 6 function guarantees that the maximum value of u is Uq. One may find that if taking 



{ujq = 2A), the above solution for the B-meson wavefunction agrees with the result in Ref. 11 1. 
Eq.(|14j) shows that the dependence on transverse and longitudinal momenta is strongly 
correlated through the combination k5^/[ci;(ci;o — ^)]- Similar transverse momentum behavior 

. Il5l and has been obtained by using the dispersion relations and 
16|, where they stated that the kr dependence of the wavefunction 



has been discussed in Re: 
the quark-hadron duality 



depends on the off-shell energy of the valence quarks, i.e. ~ ]<.\/x{l — x), where x is the 



momentum fraction carried by the valence quarks. Under such kind of /cT-dependence, 
by doing the Fourier transformation, we find that the correlation between uj and for the 
transverse part combines in a way like [z^C(j(co'o — cu)]. Thus the above results has been proved 
completely and Eqs. (j4l5p impose a strong constraint, i.e. xi'^i^'^) = ' (^o ^ ^) ' -2^], on 
the B-meson wavefunction. 

However, one may ask a question whether the ansatz (jH)) is too strong. In order to study 
the B-meson wavefunction more deeply, we assume that the B-meson wavefunction "$±{00, z^) 
can be constructed in a more general way like, 

vl/±(^, z") = <l>±(^) ■ [piz") + k{u)] ■ x[uJ ■ K - oj) ■ z^] . (16) 

Note here $±(^7) has no exact meaning of the distribution amplitude. Substituting the above 
equation into Eq.(j3]), one may find that the non-trivial solutions can only be obtained under 
the following conditions. 



p\z^) 



z'[^+{uj) - + [^+{uj) + uj^'_{uj)] = 0, (17) 



cj$_(cj)k'(cu) + ($+(cj) + lj^'_{lu))k{lu) = 0, (18) 
$+(u;) ^$_(o^) = 0. (19) 

It can be seen that these equations have the following solutions, 

p{z') = hiz' f, k{u;) = /3^-^(^o - ^)^^ (20) 
<I>+(^) = /2^'+^(^o - i^f\ <f-(^) = f2Uj''{uJo - i^y^'', (21) 

where C is arbitrary and fi{i = 1,2, 3) are undetermined parameters that have nothing to do 
with and u. Substituting these solutions into Eq.(jS)) and doing the variable transformation 

z"^ x/[uj{ujq — u)], we obtain 



with 



u^gi{x) + i8uoUJ-2co^Q)g2ix) =0, (22) 



9iix) = xhi-3x{x)-2{-6 + x)x\x)+Axi7x"ix) + 2xx^'\x))) + 
x^fi{{4C^ + 8C^ - 3x - 2Cx)x{x) + 2x((6 + AC {A + 3C) - 
x)x'{x) + 2x{{7 + 6C)x"{x) + 2xx^'\x)))), 

92{x) = xh{2x{x) + (-8 + x)x!{x) - Ax{Ax!'{x) + xx^^\x))) - 



x^/i((4C2(l + C) - (2 + C)x)x{x) + x((8 + 4C(5 + 3C) - 
x)x'{x) + 4x((4 + 3C)x"(x) + xx^''\x)))), 

where x'i^) = dx{x)/dx, x"ix) = d'^x{x) /d'^x and x^^Kx) = d'^x{x) / d'^x. The functions 
gi{i = 1,2,3) must be set to zero so as to ensure that Eq. ()22j) always be tenable with the 
variation of uj. And then we obtain that x{x) niust satisfy the following equation, 

ix'lhx"^ + fs]x"{x) + [8C/ixi+^ + 4/ix^+^ - 
4/3X] ■ x'{x) + [4CVix^ - + hxMx) = 0. (23) 

The solution of Eq. (f23|l from the condition of /s + fix'-' = should be excluded, because 
when setting /s + /ix*" = 0, one may find that it will lead to "$±{00, z"^) = 0. Under the 
condition that /s + fix'^ 7^ and by doing the transformation, x{x) — * ^(x)/(/3 + /ix*"), we 
obtain 

4xC{x) + 4^'{x) - ^{x) = 0. (24) 

Eq.(j21I) can be changed into a zeio-th normal Bessel function by doing the variable trans- 
formation, X \J —X [x < 0), and then we get x{x) = Jo{V~x)/ {f^ + fix^). 

With the help of the above results, we finally obtain the solution for the B-meson wave- 
functions under the ansatz Eq. lfTHjl . i.e. 

i>+{uj, kl) = ^d{ujQ - u)5 (ki - u{uoq - uj)) , (25) 

§_(^, ki) = '^-^t^^^^eiyj^ - u)b (ki - ^(^0 - ^)) , (26) 

where the normalization condition, / (iti;(i^k_L\E'-|-(ti;, k^) = 1, has been adopted. 



IV. PHYSICAL CONSEQUENCE AND DISCUSSION 

We have solved the equations which are derived from the equation of motion of the light 
spectator quark by using two different ansatz, i.e. Eq.® and Eq. ()16|) . Under the condition 
ijJq = 2A, the results (see Eqs. fll4pi5|l and Eqs. ()25p2(jj) ) for these two different^ansatz agree 
with the one obtained directly by taking the heavy quark limit mt cxd|1i|, i.e. these 
three prescriptions are equivalent. Eqs. ()25l26|) show that the dependence on the transverse 
and longitudinal momenta is strongly correlated through a non-factorizable combination 
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(a) (b) 
FIG. 1: The left diagram is the distributions of (f)^^{uj) and (j)^^{uj); the right diagram is the 
distribution of the ratio R^^{uj) = (uj) / (p^^ (uj) . In the left dia gram, the thicker same type 
lines are for (j)^^{uj) and the thinner ones are for (j)^^{uj), respectively. 

l^j_/[^(^o ~ ^)]- As has been pointed out in Refs. \v\ . such kind of transverse momentum 
dependence for the B-meson wavefunction has a slow-dumping with oscillatory behavior at 
large transverse distances, "^±{10, z"^) ~ cos (jz±\yJuj{ujQ — cu) — 7r/4^ / ^J\z±]. This behavior 
is quite different from that of the models with a simple Gaussian distribution, which has a 
strong dumping behavior at large transverse distances. 

Next, we make a discussion on the behavior of the obtained B-meson distribution am- 
plitudes. For latter convenience, we label the B-meson distribution amplitudes in Eq. (jl5p 
as The value of ctjn in Eq. ()15j) can be taken as the continue threshold of the light 



quarks in the B meson P. [l9j. which is usually taken to be in the interval (0.8 — 1.0)GeV^. 
Other than taking ^±(^7) directly into calculations, one usually takes the combined form of 

= = (27) 

Another typical definition for 0_b(cl') and (psi^) can be found in Ref.|2(jj, however one may 
find that the qualitative conclusions are similar. We show the distributions of (p^^iu) 



and 



with varying ujq in Fig.(^) and the distributions of the ratio R 



NBi 



(f)^^ {uj) / (j)^^ {lo) in Fig.((TjD). One may observe that the value of (f)^^{uj) is always bigger 



^ observing that A e (0.4, 0.5)Gey, ujq roughly agrees with the relation ujq = 2A. 
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FIG. 2: The left diagram is the distributions of (j)'^^{ui) and (f>'^^{u>); the right diagram is the 
distribution of the ratio R'^^{u)) = ^'^^ {uo) / (j/^^ {uo) . In the left diagram, the thicker same type 
lines are for (j)'^^{uj) and the thinner ones are for (j)^^{iAj), respectively. 

than (j)^^{uj) and the value of 0^'^(co') is negative in small u regions. 

Some other models for the B meson DAs have also been adopted in the literature. In 
Ref. Ql, the authors adopted the model, 



which are obtained under the prescription that the transverse part and the longitudinal part 
of the B wavefunction can be strictly factorizedj^ 3]- In fact, (p'^^iuj) can be derived from 
Eqs. fl4l5|) by requiring dxi'^ ■, z^) / duj = together with the ansatz Eq.(jH)). The parameter 
A can be determined through = 4dx{i^, z"^) /dz'^\;^2^Q, whose value is of order Aqc'dQi. 
Inspired by a QCD sum rule analysis, Grozin and Nuebert |17] have proposed a simple model 
for the distribution amplitudes, 

where Qq = 2A/3. Both (j)'^^{uj) and (p'^^luj) satisfy the relation Eq.((7j). It means that 
(p^^iuj), (p^^iu) and (pl^iuj) all are solutions of Eqs.(jlEl) as ^ 0. (j)^^{uj) and (pl^iu) 
have different behaviors, however they have the same asymptotic behavior that is favored 
by most of the calculations [3, 0, 



1-71 II ok J.NB,GN/ N >NB,GN/ n , ^ 

171 IIHI . , (uo) uj, (p^ [uj) ~ const, as ^ (J. 



4>^^{uj) and (p^^iuj) have similar behavior (exponential form), but the asymptotic behavior 
for 0+ is different, i.e. (p^^iui) ~ a;^, as a; ^ 0. More clearly, we show the distributions 



of and and the ratio R^^{uj) = {uj) / [uj) with varying A in Fig.©. 

From Fig. ((21), one may observe that similar to the case of and (p'^^luj), (j)'§^{uj) is 

also comparable to (f)%^{uj) in the endpoint regions. 

One may observe from Figs. ()ip2j) that the value of (j)^^''^^ {uo) is always bigger than 

— NB GS — NB GS 

(j)^ ' {u) and the value of (j)^ ' {u) is negative in small u regions, such behavior might 
lead to the total net contribution frorn_^^_f_^£(a;) be much smaller than that of 



iNB,GS, 



In literature, many authors (see Refs. 2l|, |2^, |23|) did the phenomenological analysis with a 



single distribution amplitude 0^, setting 0^ = (or strictly speaking, ignoring the contri- 

— N B GS N B GS 

butions from (ps)- However, since 0+ ' {uo) and 0_ ' {uo) have a quite different endpoint 
behavior, such difference maybe strongly enhanced by the hard scattering kernel. The re- 
sults in Refs. ^ for B ^ n transition form factor confirm this observation. Especially 



in Ref.j25|, by comparing the PQCD results in the large recoil regions with those obtained 
from the QCD light-cone sum rules and the extrapolated lattice QCD, the authors give a 
detailed analysis on the consistent calculation of the B ^ tt transition form factor in the 
whole physical region. Their results show that a better slope of the PQCD results can be 
obtained by taking both (ps and (pB into consideration. In fact, a discussion on using a single 
(pB is given in Ref.j2^, which is based on the assumption that 0_ vanishes at the both ends 
of the momentum u. There is therefore no convincing motivation for setting 0^ = and 
such an approximation may lead to unreliable results. Since for the endpoint region, the 
contribution from the hard scattering part might be big, which is the case of B ^ n and 
B — i> transition form factors, so we argue that both 0^ and 0b should be kept for a 
better understanding of the B physics. 



V. CONCLUSION 



In summary, we have solved the B meson wavefunction based on the equations derived 
from the vacuum-to-meson matrix element of the nonlocal operators in the heavy quark 
effective theory and from the equation of motion of the light spectator quark in the B 
meson. Our analysis shows that the equation of motion of the light spectator quark can 
impose a strong constraint on the B-meson wavefunctions "^±{00, z"^). For example, the 
function xi^y -2^) depends only on a single combined variable [u ■ {uq — 00) ■ z"^], if assuming 
"^±{00, z"^) = (I)±{lj)x{uj, z'^). According to our discussion, the distribution amplitude (psil^) 
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of the B meson can not be safely neglected for a better understanding of the B decays. 
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